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The  Upper  Envelope  of  Piecewise  Linear  Functions: 
Algorithms  and  Applications 

Herbert  Edelsbrunner^,    Leonidas  J.  Guibas"'    and    Micha  Sharir'* 


Abstract.  This  paper  studies  applications  of  envelopes  of  piecewise  linear 
functions  to  problems  in  computational  geometry.  Among  these  applications 
we  Hnd  problems  involving  hidden  line/surface  elimination,  motion  plan- 
ning, transversals  of  polytopes,  and  a  new  type  of  Voronoi  diagram  for  clus- 
ters of  points.  All  results  are  either  combinatorial  or  computational  in 
nature.  They  are  based  on  the  combinatorial  analysis  in  two  companion 
papers  (Pach  and  Sharir  ;PSj  and  Edelsbrunner  E2|)  and  a  divide-and- 
conquer  algorithm  for  computing  envelopes  described  in  this  paper. 
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1.    Introduction 

This  paper  continues  the  study,  initiated  in  [PSl  and  in  (E2],  of  envelopes  of  piecewise 
linear  functions  in  two  or  more  variables.  The  previous  papers  have  established  tight 
lower  and  upper  bounds  on  the  combinatorial  complexity  of  such  envelopes.  In  this 
paper  we  provide  efficient  algorithms  for  calculating  envelopes  of  this  kind,  discuss 
several  extensions  and  special  cases  of  the  previous  combinatorial  bounds,  and  give  a 
variety  of  applications  of  these  results  to  many  problems  in  combinatorial  and  compu- 
tational geometry. 

Let  us  first  review  the  results  of  PS!  and  'E2|.  Let  Si,52,...,s„  be  n  d-simplices 
in  (</+l)-dimensional  space,  none  of  which  is  vertical  (that  is,  parallel  to  the  (<^  +  l)'' 
coordinate  axis).  We  can  thus  view  each  5,  as  the  graph  of  a  partially  defined  linear 
function  x^^i=  f,{xi,X2,...,JC^),  whose  domain  of  definition  is  a  <f -simplex,  namely  the 
orthogonal  projection  of  s,  onto  the  hyperplane  x^^i=0.  The  upper  envelope,  M,  of 
the  given  simplices  is  the  pointwise  maximum  of  these  functions,  that  is, 

.U(ii,X2,...,Jrf)  =    max   |/,(ji,X2,...,Xj)|, 

where  each  /,  is  assumed  to  be  — x  outside  its  domain  of  definition.  The  tower 
envelope  of  the  simplices  is  defined  in  a  symmetric  fashion. 

One  can  associate  with  the  envelope  .V/  a  polyhedral  cell  complex,  .V/*,  in  <f -space 
such  that  over  each  cell  of  M*  the  envelope  M  is  attained  by  a  fixed  function  /,. 
Intuitively,  this  is  the  orthogonal  projection  of  the  graph  of  M  onto  x^^i=0  (see  Fig- 
ure 2.1  which  shows  the  projection  of  four  triangles  in  d+l=3  dimensions).  The 
combinatorial  complexity  of  M  is  the  complexity  of  this  complex,  that  is,  the  total 
number  of  faces  (of  any  dimension)  composing  it.  In  general,  the  projection  of  the 
faces  of  M  does  not  yield  a  convex  decomposition  of  «f-space.  However,  one  can 
obtain  a  refined  convex  decomposition  by  superimposing  M*  onto  the  arrangement  of 
the  (((/  — l)-dimensional)  hyperplanes  containing  the  (d— l)-faces  of  the  given  sim- 
plices (see  Figure  2.2  which  shows  M*  for  a  set  of  three  triangles  in  d+l=3  dimen- 
sions).   We  will  measure  the  complexity  of  M  in  terms  of  this  refined  decomposition. 

If  instead  of  simplices  we  have  a  collection  of  arbitrary  piecewise-linear  functions 
of  d  variables,  we  can  decompose  the  graph  of  each  of  them  into  a  collection  of  sim- 
plices, and  then  obtain  the  upper  envelope  of  the  given  functions  as  the  upper 
envelope  of  these  simplices. 

The  mentioned  two  previous  papers  analyze  the  combinatorial  complexity  of  such 
envelopes  in  d  +  l  dimensions.  They  show  that  it  is  0(n  a{n)),  where  a{n)  is  the 
extremely  slowly  growing  inverse  of  Ackermann's  function.  Moreover,  this  bound  is 
tight  in  the  worst  case.  For  d+l  =2,  we  face  the  special  case  of  the  envelope  of  n  line 
segments  in  the  plane.  This  case  has  been  studied  in  [HS|,  [WS],  [Si,  where  it  is 
shown  that  the  complexity  of  the  envelope  (in  this  case  the  number  of  subsegments 


composing  it)  is  0(na{n)),  and  that  this  bound  is  tight  in  the  worst  case.  The  proofs 
are  based  on  reformulating  the  problem  in  terms  of  Davenport  Schinzel  sequences  of 
order  3. 

As  a  matter  of  fact,  in  the  one-dimensional  case,  the  theory  of  Davenport- 
Schinzel  sequences  yields  tight  almost  linear  upper  bounds  on  the  complexity  of  the 
envelope  of  any  collection  of  (partially  defined)  continuous  univariate  functions,  pro- 
vided that  each  pair  of  them  intersect  in  at  most  a  fixed  number  of  points.  In  con- 
trast, for  collections  of  n  functions  of  d  >2  variables  (which  satisfy  appropriate  condi- 
tions on  the  pattern  of  their  intersections),  no  comparable  tight  upper  bounds  on  the 
complexity  of  their  envelopes  is  known  in  general  as  yet,  and  the  only  general  bound 
available  so  far  is  the  trivial  one,  0(n  "').  We  refer  to  [SSj  for  a  few  improved 
results  for  certain  classes  of  bivariate  functions. 

The  proofs  of  the  mentioned  bounds  for  piecewise  linear  functions  use  induction 
on  the  dimension  d.  The  proof  in  'PSi  is  based  on  a  divide-and-conquer  analysis. 
One  partitions  the  collection  of  n  simplices  into  two  subcollections  of  roughly  n/2 
simplices  each,  recursively  obtains  the  envelope  of  each  subcollection.  and  then 
obtains  the  overall  envelope  by  taking  the  pointwise  maximum  of  the  two  suben- 
velopes.  Using  arguments  based  on  arrangements  of  hyperplanes,  convexity,  and  on 
the  complexity  of  the  envelope  restricted  to  certain  lower  dimensional  spaces,  one 
shows  that  the  number  of  additional  facets  created  by  superimposing  the  two  suben- 
velopes  is  bounded  by  0(n  a{n)).  This  leads  to  a  recurrence  relation  whose  solution 
yields  the  desired  bound  on  the  number  of  facets  (that  is,  <i-dimensional  faces)  in  the 
envelope.  Using  Euler's  formula  for  planar  maps  this  yields  similar  bounds  on  the 
total  number  of  faces  of  M  if  d+l—3.  For  higher  dimensions  this  proof  has  been 
extended  in  E2j  who  obtains  the  same  bound  for  the  overall  complexity  of  the 
envelope  using  a  different  divide  step. 

In  Section  2,  we  follow  the  outline  of  the  proof  in  [PSj  to  obtain  a  divide-and- 
conquer  algorithm  for  calculating  the  envelope  of  triangles  in  three  dimensions.  The 
amount  of  time  needed  is  at  most  proportional  to  the  maximum  combinatorial  com- 
plexity of  the  envelope,  namely  0(n  oi{n)).  Because  of  the  lack  of  convex  hull  algo- 
rithms in  four  and  higher  dimensions  that  run  in  time  proportional  to  their  output, 
we  have  not  succeeded  in  generalizing  this  algorithm  so  that  it  runs  in  time 
0(r  ot(n))  if(/+l>4.  We  also  discuss  several  extensions  and  special  cases  of 
envelopes  of  piecewise  linear  functions  which  are  needed  for  the  applications  that  we 
study. 

The  major  part  of  the  paper  is  devoted  to  applications  of  the  combinatorial  and 
computational  bounds  stated  above.  These  applications  include  the  hidden 
line/surface  removal  problem  which  is  discussed  in  Section  3.  We  give  algorithmic 
results  that  match  and  generalize  previous  results  of  [Dv|,  [Mj  for  d+l=3  dimensions. 
We  also  obtain  algorithms  for  related  problems,  such  as  constructing  an  image  of  a 
solid    defined    in    Constructive    Solid    Geometry,    and    obtaining    views   of   a    three- 


dimensional  projection  of  tetrahedra  in  four  dimensions. 

Section  4  considers  translational  motion  planning  for  polyhedra  in  three  dimen- 
sions. Here,  we  calculate  the  space  of  free  placements  of  a  given  polyhedron,  B, 
which  is  free  to  translate  amidst  a  collection  of  polyhedral  obstacles.  We  also  discuss 
special  cases  of  the  problem,  such  as  where  the  obstacles  form  a  polyhedral  terrain 
(that  is,  a  piecewise  linear  surface  meeting  each  vertical  line  at  exactly  one  point),  and 
the  case  where  B  is  allowed  to  translate  from  its  current  position  along  a  straight 
half-line  only.  The  latter  special  case  extends  work  on  movable  separability  of  sets 
reported  in  [Tt]. 

The  problem  of  stabbing  line  segments  and  polytopes  is  investigated  in  Section  5. 
Using  a  standard  duality  transformation,  we  show  that  stabbing  hyperplanes  can  be 
represented  by  points  lying  between  the  upper  envelope  of  one  collection  of  simplices 
and  the  lower  envelope  of  another  such  collection.  Our  results  extend  previous  work 
on  this  problem  in  two  dimensions  (see  [E3|). 

Voronoi  diagrams  of  point  clusters  are  considered  in  Section  6.  For  each  cluster 
(that  is,  finite  set  of  points  in  d  dimensions)  its  distance  from  a  point  is  the  maximum 
Euclidean  distance  from  this  point  to  any  point  in  the  cluster.  The  Voronoi  diagram 
of  a  collection  of  clusters  is  then  defined  as  the  decomposition  of  space  into  maximal 
domains  so  that  for  each  domain  there  is  a  unique  nearest  cluster  for  all  points  in  this 
domain.  By  transforming  the  problem  to  ((f -|-l)-dimensional  space,  we  can  reformu- 
late it  in  terms  of  envelopes  of  certain  piecewise  linear  functions. 

We  do  not  regard  this  list  of  applications  as  exhaustive,  and  we  expect  many 
more  applications  to  be  found.  For  example,  Tamir  [Tr]  has  recently  discovered 
applications  of  our  results  to  p-center  and  obnoxious  p-center  problems  for  certain 
trees  and  graphs. 


2.    Computing  an  Envelope 

This  section  presents  an  algorithm  for  constructing  the  upper  envelope  of  a  set  of  n 
triangles  in  three  dimensions.  The  algorithm  follows  the  outline  of  the  proof  in  PS| 
that  shows  that  the  combinatorial  complexity  of  this  envelope  is  0(n^a(n)).  At 
several  points  we  have  to  introduce  intricate  algorithmic  tools  in  order  to  get  a 
worst-case  optimal  algorithm.  For  some  of  these  tools  the  complexity  goes  up  more 
than  desired  when  we  generalize  them  to  four  and  higher  dimensions.  This  explains 
why  we  do  not  have  an  optimal  (or  even  near-optimal)  method  for  computing 
envelopes  in  four  or  higher  dimensions  yet.  After  presenting  and  analyzing  the  algo- 
rithm, we  study  a  few  extensions  of  envelope  constructions.  These  will  lead  to  several 
computational  and  combinatorial  results  used  in  later  sections  of  this  paper. 

We  next  present  the  algorithm  that  constructs  the  upper  envelope  of  a  set,  S,  of 
n  triangles  in  three  dimensions.  Whenever  convenient  in  the  discussion  we  will  make 
implicit  assumptions  about  the  triangles  being  in  general  position.    The  main  reason  is 


that  we  hope  to  get  the  point  across  if  we  leave  out  tedious  complications.    We  see  the 
general   method,  called   the   "simulation  of  simplicity",   described   in    [EM}    (see  also 
El  ),  as  a  justification  of  this  sloppy  attitude. 

First,  we  need  a  few  definitions.  We  write  M{S]  for  the  upper  envelope  of  5, 
and  M*{S)  for  the  subdivision  obtained  by  projecting  the  faces  of  M{S)  vertically 
onto  the  plane  ^3  =  0  (see  Figure  2.1  which  is  borrowed  from  PS  ).  In  general  the 
regions  of  .U*(S)  are  not  convex.  To  make  them  convex  we  refine  M*{S)  by  project- 
ing   all    triangles    vertically   onto    X3=0    and    extending    the    3|5|    triangle    edges    to 

unbounded  lines.  The  arrangement^  defined  by  these  lines  is  denoted  by  A  (5),  and 
A/(S)  denotes  the  refined  subdivision  that  we  get  by  superimposing  M*{S)  and  A (5). 
See  Figure  2.2  taken  from  [PS|;  it  shows  the  projection  of  three  triangles  and_the 
extension  of  their  edges  yielding  an  arrangement  of  nine  lines.  Jn  order  to  make  .V/(5) 
a  viable  representation  of  \f{S)  we  associate  each  region  of  -V/(S)  with  a  gointer  to 
the  triangle  that  assumes  the  maximum  height  over  this  region.  Since  \t{S)  is  a 
refinement  of  M^{S)  this  triangle  in  unique  in  any  case.  .V/*(S)  and  also  \'f(S)  can  be 
obtained  from  -V/(5)  by  merging  adjacent  regions  above  which  the  maximum  height  is 
assumed  by  the  same  triangle.  By  construction,  the  envelope  vertically  above  a  region 
of  A{S)  is  a  concave  function. 

Algorithm  1    (constructs  M(5_)_as  a  representation  of  the  upper  envelope  of  5.) 
if  |5|<1  then  Construct  M{S)  directly    else 

Step  1.    Partition  5  into  sets  5j  and  S2  of  sizes  l.|'S|/2j  and  [|5|/2l. 


Figure  2.1. 


The  arrangement  defined  by  a  finite  set  of  lines  in  the  plane  is  the  subdivision  of  the  plane  that 
we  get  by  drawing  the  lines.  It  consists  of  vertices  (points  where  lines  intersect),  edges  (pieces  of  lines 
that  connect  vertices),  and  regions  (connected  components  of  the  plane  reduced  by  all  lines). 


Figure  2.2. 

Step  2.    Construct  A/{S_J  and  .V/(52)  recursively. 

Step  3.  Superimpose  .V/(Si)  and  .4(52)  ^^^'  symmetrically,  superim- 
pose XfiS-y)  and  .4(Si).  We  denote  the  thus  created  subdivisions  by 
M{Si)  and  \I{So). 

Step    4.     Construct    A_[S)    which    is    i4(Si)    and    .4(52)   superimposed.  . 
Thus,  M{Si)  and  M(S2)  are  refmements  of  A  (5). 

Step  5.  For  each  region  r  of  A (5)  and  for  i  =  1,2  construct  set  S,  rC5, 
that  containj_all  triangles  of  5,  assuming  the  maximum  height  over 
a  region  of  .V/(5,)  contained  in  r. 

Step  6.  For  each  region  r  of  .4(5)  construct  .V/(5)  restricted  to  r  by 
intersecting  the  half-spaces  bounded_from  below  by  the  planes  that 
contain  the  triangles  in  5i  rU52  r.  M(5)  restricted  to  r  is  the  verti- 
cal projection  of  the  boundary  facets  of  this  convex  polyhedron 
clipped  at  the  boundary  of  r. 
endif. 

Let  us  now  discuss  the  various  steps  of  .\lgorithm  I  in  further  detail.  At  the 
same  time  we  analyze  the  time-complexity  of  each  individual  step  which  will  then 
lead  to  the  recurrence  relation 


T{n)=2T(f)+0{n'a(n)), 

for  the  amount  of  time,  T{n),  the  algorithm  takes.    This  recurrence  relation  solves  to 
0{n^a{n))  (as  in  [PS|). 

The  first  non-trivial  step  of  Algorithm  I  is  Step  3  which  superimposes  \[(S,)  and 
'4(53_,),  1=1,2.  The  combinatorial  analysis  in  !PSl  guarantees  that  the  number  of 
faces  in  the  resulting  subdivision,  .V/(5,),  is  0{n-a{n)),  n  =  |5i|-t-|52|.    Since  all  regions 


in  both  subdivisions,  A/(5,)  and  .4(53_,),  are  convex  we  can  use  the  superimposition 
algorithm  of  Guibas  and  Seidel  GS;  which  takes  time  linear  in  the  size  of  the  output. 
Thus,  Step  3  takes  0{n^a{n))  time. 

Step  4  superimposes  two  arrangements  which  can  be  done  in  0(n'')  time. 
n=|5i|+|52|,  using  again  the  same  superimposition  algorithm.  Alternately,  we  could 
construct  the  resultant  arrangement,  A  (5),  from  scratch  which  also  takes  only  0{n'] 
time  (see  [El;). 

To  understand  Step  5  it  is  important  to  recall  that  subdivision  .V/(5,),  for  i'  =  l,2, 
as  constructed  in  Step  3  is  a  refmement  of  arrangement  A{S)  which  is  constructed  in 
Step  4.  Thus,  e^ch  region  of  A{S)  contains  several  regions  of  >/(5,).  Recall  also  that 
each  region  of  M{S,)  records  the  triangle  that  assumes  the  upper  envelope  above  this 
region.  The  goal  of  Step  5  is  to  colject,  for  each  region  r  of  A{S),  the  triangles  associ- 
ated with  regions  in  -V/fSJ  and  -V/fSj)  contained  in  r.  This  can  be  done  by  visiting 
all  regions  of  M{S,)  using  a  graph  traversal  algorithm  that  first  exhausts  all  regions 
contained  in  a  common  region  of  the  arrangement  and  then  goes  to  an  adjacent 
arrangement  region.  This  is  a  straightforward  application  of  depth-first  search_(see 
,Tn  )  which  takes  time  proportional  to  the  number  of  regions  and  edges  of  A/(5, ). 
Thus,  Step  5  also  takes  time  0(n  a{n)). 

Finally,  we  consider  Step  6  which  is  the  most  subtle  part  of  Algorithm  I.  For 
each  region  r  of  A{S)  it  constructs  the  intersection  of  the  half-spaces  that  are 
bounded  from  below  by  the  planes  that  contain  the  triangles  in  Sj  rU-?,  r-  ^f  ''  is  the 
J  region  of  .4(5)  and  f",  =  |5,  J+ISj  ^ |,  then  this  can  be  done  in  O(m,logm,)  time 
(see  e.g.  'PrS|).  Unfortunately,  we  only  know  that  the  sum  of  the  m,  is  0{n'^a{n)) 
which  does  not  imply  anything  better  than  that  the  sum  of  the  m,  logjm,  is 
0(n'^a{n)\ogn).  Thus,  to  achieve  0(n^a(n))  running  time  we  have  to  intersect  the 
half-spaces  in  a  somewhat  more  intelligent  manner.  To  describe  such  an  alternate 
algorithm  define  the  slope  of  a  triangle  as  the  slope  of  the  line  of  intersection  between 
the  plane  spanned  by  the  Xj-  and  the  Xj-axis  and  the  plane  that  contains  the  trian- 
gle, [f  the  slopes  of  all  triangles  in  S,^  were  smaller  than  the  slopes  of  all  triangles  in 
S2r  we  could  construct  the  intersection  of  the  half-spaces  by  merging  the  recursively 
constructed  polyhedra  for  Sj  ,  and  5,  ^  in  linear  time  (see  PrSl).  Unfortunately 
again,  the  recursion  is  such  thai  the  polyhedra  constructed  are  not  exactly  those  for 
the  regions  of  .4(5)  but  for  the  regions  of  A{Si)  and  A (52).  For  example,  let  r,  be 
the  region  of  -4(5,)  that  contains  region  r,  for  i'  =  l,2.  It  is  not  advisable  to  use  the 
polyhedron  representing  the  subenvelope  above  r,  as  a  substitute  for  r's  polyhedron 
since  it  might  have  many  facets  that  belong  to  half-spaces  redundant  above  r.  The 
combinatorial  bound  on  the  sum  of  the  m,  does  not  generalize  to  these  larger 
numbers;  thus  it  is  crucial  not  to  be  generous  at  this  point. 

The  way  out  of  this  dilemma  is  to  remember  that  the  subdivision  of  r,  in  \f{S,) 
can   be  viewed   as   a   representation  of  the   polyhedron  of  r, .     The  separation  of  the 


triangle  slopes  implies  that  the  intersection  of  the  boundaries  of  the  two  polyhedra, 
the  ones  of  r,  and  rj  restricted  to  the  area  above  region  r  in  A{S),  is  a  connected  and 
piecewise  linear  curve.  Figure  2.3  displays  r^  rj,  r,  and  the  curve  without  showing 
the  decompositions  of  the  regions.  Keep  in  mind,  however,  that  this  curve  can  merge 
into  the  boundary  of  r  and  leave  it  again  an  arbitrary  number  of  times.  Because  of 
the  slope  condition,  this  curve  has  the  property  that  it  intersects  any  plane  normal  to 
the  X2-axis  in  a  single  point.  Using  standard  methods  for  merging  two  subdivisions 
along  a  monotone  curve  (see  [PrS]  and  papers  on  merging  Voronoi  diagrams  referred 
to  in  iPrSj),  the  total  amount  of  effort  i|Jinear  in  the  number  of  edges  of  r  plus  the 
number  of  regions  that  subdivide  r  in  M(Si)  and  ^(Sj).  In  order  to  make  this  all 
work  we  have  to  provide  the  appropriate  subdivisions  of  the  regions  of  A(S).  But 
these  are  provided  by  the  superimposition  of  iV/(S,)  and  •A(53_,)  which  decomposes 
the  subdivision  of  r,  into  smaller  pieces  coinciding  with  regions  in  A(S).  Thus,  this 
superimposition  in  Step  3  implicitly  constructs  the  proper  polyhedra  (or  suitable 
representations  by  subdivisions  thereof)  which  can  then  be  merged  in  linear  time  each. 
The  only  unresolved  problem  is  now  how  we  can  guarantee  that  the  slopes  of  the 
triangles  in  Sj  ^  are  smaller  than  the  slopes  of  the  triangles  in  52^-  But  this  can  be 
achieved  if  the  initial  partitioning  step  of  Algorithm  1  constructs  S^  and  52  intelli- 
gently rather  than  arbitrarily.  Just  take  the  [|5|/2j  triangles  with  smallest  slopes, 
call  this  set  Sp  and  define  52  =  S— S,. 

We  thus  have  an  optimal  algorithm  for  constructing  the  upper  envelope  of  n  tri- 
angles in  three  dimensions.  This  is  .\lgorithm  1  with  two  changes.  First,  the  parti- 
tion of  the  set  of  triangles  takes  into  account  the  slopes  of  the  triangles.  Second,  Step 
5  is  now  superfluous  and  can  be  removed.  This  implies  the  main  result  of  this  sec- 
tion. 

Theorem  2.1.  The  upper  envelope  of  a  set  of  n  triangles  in  three  dimensions 
can  be  constructed  in  0(ra  oi{n))  time  and  storage.  This  is  optimal  in  the 
worst  case. 


-x^ 


Figure  2.3. 


We  remark  that  Algorithm  1  can  be  modified  so  that  it  constructs  the  upper 
envelope  of  n  line  segments  in  two  dimensions  in  0(na(n)logn)  time  and  0{na{n)) 
storage.  The  amount  of  storage  is  optimal  since  the  envelope  can  consist  of  Q{na(n)) 
edges;  whether  or  not  the  time  bound  is  optimal  is  still  an  open  problem.  There  is  no 
difficulty  in  generalizing  Algorithm  1  to  four  and  higher  dimensions,  however,  it  is 
still  an  open  problem  whether  or  not  this  can  be  done  such  that  the  running  time  is 
0(n  oc{n))  for  n  (f-simplices  in  d  +  i  dimensions.  This  would  then  be  optimal  since 
the  combinatorial  complexity  of  the  envelope  is  Q{n  a(n))  in  the  worst  case.  The 
main  obstacle  in  obtaining  this  result  is  Step  6  which  intersects  half-spaces. 
Currently  there  is  no  algorithm  available  that  takes  less  than  Q(m  )  time,  where  m  is 
the  number  of  half-spaces,  no  matter  how  many  or  few  faces  the  resulting  polyhedron 
has.  The  combinatorial  results  in  PSj  and  E2|  only  bound  the  sum  of  the  m,  (where 
m,  is  the  number  of  non-redundant  half-spaces  above  the  i  cell  of  the  d~ 
dimensional  arrangement)  and  not  the  sum  of  the  m,  .  Indeed,  there  are  cases  where 
the  sum  of  the  m,  is  Q(n  ~  )  and  thus  contradict  the  desired  0{n  a{n))  upper 
bound.  An  approach  that  might  be  worth  pursuing  is  to  design  an  algorithm  that  fol- 
lows the  outline  of  the  divide-and-conquer  proof  in  [E2'.  The  main  difference 
between  Algorithm  1  and  such  a  hypothetical  algorithm  would  be  that  the  latter 
recurs  for  a  constant  number  of  non-disjoint  subsets  of  (i-simplices  rather  than  for 
two  disjoint  subsets. 

The  remainder  of  this  section  studies  three  e.xtensions  of  the  envelope  problem 
which  have  algorithmic  as  well  as  combinatorial  applications  later  in  the  paper.  The 
first  extension  considers  the  region  of  points  that  lie  above  the  upper  envelope  of  a 
finite  set  5  of  <f-simplices  and  below  the  lower  envelope  of  another  finite  set  T  of  d- 
simplices  in  d  +  l  dimensions.  From  the  combinatorial  results  in  [PS  and  E2  we 
know  that  both  envelopes  have  complexity  0(n  a{n)),  with  a=|5|  +  |r|,  and  by 
Theorem  2.1  we  can  construct  both  envelopes  in  0{n^a{n))  time  if  d  +  l  =3.  But  how 
can  we  be  sure  that  the  intersection  of  the  two  envelopes  does  not  exceed  these  com- 
plexity bounds?  In  fact,  it  does  not.  One  way  to  see  this  is  to  go  through  the  proofs 
of  the  0(n  ci{n))  upper  bounds  for  upper  envelopes  and  to  make  one  crucial  change: 
rather  than  constructing  the  two  envelopes  above  a  cell  of  the  (/-dimensional  arrange- 
ment separately,  we  construct  the  region  of  points  between  the  two  envelopes.  Res- 
tricted to  the  area  above  a  cell  r  of  the  arrangement,  this  region  is  the  intersection  of 
half-spaces  and  thus  convex.  It  follows  that  its  complexity  is  bounded  by  the  sum  of 
the  complexities  of  the  two  regions  between  the  two  corresponding  upper  and  lower 
subenvelopes.  This  is  all  we  need  to  get  the  desired  combinatorial  result;  also  Algo- 
rithm I  still  works  nicely  in  this  extended  case  if  rf-i-l  =3. 

The  second  extension  considers  the  special  case  where  each  (i-simplex  in  d  +  l 
dimensions  is  a  half-hyperplane,  that  is,  a  portion  of  a  hyperplane  h  restricted  to  one 
side  of  a  ((f  — l)-flat  in  h.  In  d  +  l=2  dimensions  a  half-hyperplane  is  a  half-line.  For 
n  such  half-lines  it  is  not  difficult  to*show  that  the  number  of  edges  in  the  upper 
envelope  is  at  most  2n  (see   E3  ).    This  two-dimensional  result  can  now  be  used  as  the 


base  case  of  the  inductive  analysis  of  upper  envelopes  in  higher  dimensions.  Recall 
that  the  only  reason  for  the  a{n)  factor  in  the  complexity  of  general  upper  envelopes 
is  that  the  base  case  considers  line  segments,  and  the  upper  envelope  of  line  segments 
has  worst-case  complexity  Q{na{n)).  The  reduction  for  half-hyperplanes  in  d  +  l 
dimensions  leads  to  a  linear  number  of  sets  of  half-hyperplanes  in  d  dimensions  whose 
upper  envelopes  have  complexity  0(n  ~^)  by  inductive  assumption.  This  leads  to  an 
upper  bound  of  0(n'')  for  the  complexity  of  n  half-hyperplanes  in  d+l  dimensions. 
Note  that  this  bound  also  holds  for  the  combinatorial  complexity  of  the  region  of 
points  above  the  upper  envelope  of  one  set  of  half-hyperplanes  and  below  the  lower 
envelope  of  another  set  of  half-hyperplanes.  Furthermore,  Algorithm  1  takes  only 
0(a")  time  if  its  input  consists  of  n  half-planes  in  three  dimensions.  This  is  because 
the  only  step  where  the  a{n)  factor  sneaks  in  (Step  3)  has  now  complexity  O(n^). 

It  is  interesting  to  note  the  similarity  between  the  upper  envelope  of  a  set  of 
half-hyperplanes  and  the  so-called  zone  of  a  hyperplane  in  an  arrangement  of  hyper- 
planes  in  d  +  l  dimensions  (see  El.  chapter  5i).  In  both  cases,  the  combinatorial  com- 
plexity is  0{n'^)  but  the  known  proofs  of  those  two  results  are  very  different. 

Finally,  we  consider  the  case  where  the  n  <f-simplices  in  d  +  l  dimensions  are 
pairwise  disjoint  (assuming  they  are  relatively  open).  In  this  case,  the  maximum 
height  above  a  cell  of  the  rf-dimensional  arrangement  is  assumed  by  only  one  d- 
simplex.  The  combinatorial  complexity  of  this  arrangement  is  0(n  )  which  implies 
the  same  upper  bound  for  the  envelope.  The  more  dramatic  effect  of  the  non- 
intersection  assumption  is  that  it  simplifies  Algorithm  1  significantly  and  thus  allows 
us  to  generalize  it  to  higher  dimensions_without  loss  of  worst-case  optimality.  Steps  3 
and  4  are  now  the  same  since  M{Sy)  =  A(Si),  .V/fSj)  =  A(S2),  and  therefore 
M{Si)  =  M(S2)  =  A{S).  In  Step  5,  set  S, ^  is  the  singleton  set  that  contains  the 
highest  simplex  in  5,  above  the  region  in  A{S,)  containing  region  r  in  A{S).  Thus, 
Step  6  simplifies  to  a  comparison  between  the  simplex  in  S^^  and  the  one  in  52  r.  for 
every  r.  The  most  expensive  step  of  this  algorithm  is  now  the  superimposition  of 
arrangements  A(5i)  and  A(S2)  which  can  be  done  in  quadratic  time  if  d  +  l=3.  In 
arbitrary  d+l  dimensions,  this  operation  takes  0{n  )  time,  for  n  =  |5il-|-|52|  (see 
El  ).  Thus,  we  have  now  an  algorithm  that  runs  in  arbitrary  dimensions  and  takes 
0{n'^)  time. 

A  similar  effect  (namely  that  the  combinatorial  complexity  of  the  envelope  is 
0{n'^)  rather  than  0{n'^a(n)))  can  be  observed  when  the  n  ci-simplices  in  d+l 
dimensions  intersect  in  a  certain  restrictive  manner.  For  example,  if  intersections 
occur  only  on  the  highest  level  of  the  recursion  (talking  in  terms  of  the  divide-and- 
conquer  algorithm)  then  this  is  true.  This  proves  that  the  upper  envelope  has  com- 
plexity 0(ri'^)  if  the  set  of  <f-simplices  is  the  union  of  two  sets  with  the  property  that 
any  two  rf-simplices  in  the  same  set  are  pairwise  disjoint.  Unfortunately  though,  the 
computational  complexity  in  this  case  might  deteriorate  to  0(n  logn),  d+l=3,  since 
the  first  divide  step  is  now  predetermined  and  cannot  separate  the  triangles  by  their 
slopes. 
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We  summarize  these  results. 

Theorem  2.2.  Let  5  be  a  set  of  n  (/-simplices  in  d  +  l  >3  dimensions  and  let  T 
be  another  such  set  whose  cardinality  is  at  most  n. 

(i)  The  region  of  points  above  the  upper  envelope  of  5  and  below  the 
lower  envelope  of  T  has  combinatorial  complexity  0{n  Oi(n)).  In  d-\-l=3 
dimensions  it  can  be  constructed  in  0(n  a{n))  time  and  storage. 

(ii)  If  all  rf-simplices  in  S  are  half-hyperplanes,  then  the  combinatorial 
complexity  of  the  upper  envelope  of  S  is  0{n  ),  and  it  can  be  constructed  in 
0(n  )  time  if  d  +  l=3.  If  also  all  (i-simplices  in  T  are  half-hyperplanes, 
then  the  same  complexity  bounds  hold  for  the  region  of  points  above  the 
upper  envelope  of  S  and  below  the  lower  envelope  of  T. 

(iii)  If  the  d-simplices  in  5  are  pairwise  disjoint,  then  the  combinatorial 
complexity  of  the  upper  envelope  is  0(n  )  and  it  can  be  constructed  in 
0{n'^)  time. 

(iv)  If  5  is  the  union  of  two  disjoint  sets  with  the  property  that  any 
two  (/-simplices  in  the  same  set  are  disjoint,  then  the  combinatorial  com- 
plexity of  the  upper  envelope  is  0{n   ). 

The  remainder  of  this  paper  considers  applications  of  Theorems  2.1  and  2.2  to 
several  problems  in  computational  and  combinatorial  geometry. 


3.    Hidden  Line  and  Surface  Removal 

Imagine  that  we  take  a  picture  of  a  three-dimensional  scene  from  a  point  at  infinity. 
To  compute  what  this  picture  looks  like  -  assuming  opaque  objects  -  is  commonly 
known  as  the  hidden  line/surface  removal  problem.  Because  of  the  importance  of  this 
problem  for  practical  applications  there  are  many  algorithms  in  the  literature  that 
were  suggested  for  the  problem  (see  e.g.  ^SSSj  for  a  classification  of  several  such  algo- 
rithms). It  is  usually  assumed  that  the  objects  in  the  scene  are  determined  by  their 
piecewise  linear  boundaries  (they  are  polytopes)  and  that  they  do  not  intersect.  We 
would  like  to  mention  that  there  is  no  essential  difference  between  the  view  from  a 
point  at  infinity  (a  parallel  view)  and  the  view  from  a  finite  point  (a  perspective  view). 
Indeed,  for  every  plane  through  a  finite  view  point  there  is  a  projective  transforma- 
tion that  maps  the  plane  to  the  infinite  plane  and  therefore  the  view  point  to  infinity. 
The  polytopes  are  mapped  to  polyhedra  such  that  the  parallel  view  from  the  new 
point  is  equal  to  the  perspective  view  from  the  old  point  on  one  side  of  the  plane. 
Note,  however,  that  this  transformation  moves  points  to  infinity  and  lets  them  come 
back  on  the  other  side  of  space.  In  other  words,  polytopes  gradually  disappear  on  one 
side  and,  at  the  same  time,  come  back  on  the  other  side  of  space.  Thus,  in  order  to 
get  a  valid  picture  we  need  to  make  sure  that  the  view  point  is  shielded  from  the 
polytopes  that  come  into  the  picture  by  traveling  through  infinity.    Such  a  shielding 
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mechanism  is  provided  if  we  map  the  original  infinite  plane  to  a  finite  plane,  using 
the  same  transformation,  and  use  this  plane  as  a  background  screen  when  we  take  the 
picture. 

In  this  section  we  adopt  a  generalized  definition  of  the  hidden  line /surface  remo- 
val problem  which  is  neither  restricted  to  three  dimensions  nor  to  non-intersecting 
objects.  We  first  discuss  the  more  standard  case  of  non-intersecting  objects  and  later 
we  extend  the  analysis  to  handle  intersecting  objects.  We  also  give  some  applications 
for  this  extension. 

In  d+l  dimensions,  the  objects  in  the  scene  are  modeled  by  a  collection  of  d- 
dimensional  simplices;  for  convenience  we  assume  that  they  are  relatively  open.  The 
problem  is  now  the  same  as  computing  the  upper  envelope  of  the  (i -simplices  assum- 
ing that  the  viewpoint  is  in  the  direction  of  the  positive  {d+l)'  coordinate  axis.  We 
can  thus  use  the  algorithms  of  Section  2  to  solve  the  hidden  line/surface  removal 
problem.  Since  we  presently  assume  that  no  two  d-simplices  intersect  (although  their 
relative  boundaries  might  intersect  which  is  important  to  allow  if  we  model  a 
polytope  by  rf-simplices),  we  can  use  Theorem  2.2  (iii)  to  get  the  following  result. 

Theorem  3.1.  Let  S  be  a  set  of  n  pairwise  non-intersecting  relatively  open  d- 
simplices  in  d+l  dimensions.  The  combinatorial  complexity  of  a  view  is 
0(n'')  which  is  best  possible  in  the  worst  case.  Furthermore,  it  can  be  con- 
structed in  0(n''-|-nlogn)  time. 

In  three  dimensions,  that  is,  if  <f+l=3,  the  same  time  complexity  was  previously 
obtained  by  'Dv|  and  'Mj  who  use  known  algorithms  for  constructing  arrangements  in 
the  plane. 

Note  that  the  0(n  )  bound  for  the  combinatorial  complexity  holds  even  if  we 
make  the  cf-simplices  translucent.  Rather  than  computing  only  the  topmost  d- 
simplex  above  a  given  point  we  determine  the  topmost  i,  for  some  constant  ^,  and  the 
"color"  at  this  point  is  a  function  of  all  £  simplices. 

Consider  next  generalizations  of  the  hidden  line/surface  removal  problem  that 
arise  when  the  given  <f-simplices  are  allowed  to  intersect.  In  this  case,  the  problem  is 
exactly  the  envelope  problem  studied  in  Section  2.  As  an  example  where  this  exten- 
sion is  needed  we  mention  an  operation  that  is  useful  in  visualizing  a  four- 
dimensional  scene  given  by  n  pairwise  disjoint  relatively  open  tetrahedra.  Project 
these  tetrahedra  onto  three  dimensions  and  compute  various  views  of  this  three- 
dimensional  scene.  Since  we  lose  one  dimension  when  we  go  from  four  to  three 
dimensions,  the  tetrahedra  in  three  dimensions  will,  in  general,  intersect.  An  alter- 
nate interpretation  of  this  operation  is  that  we  compute  views  of  the  four-dimensional 
scene  by  moving  vertical  lines  in  a  given  direction  until  they  hit  an  object.  The 
"view"  shows  the  first  object  hit  by  any  such  line.  Note  that  this  visualization  of  the 
four-dimensional  scene  is  different  from  a  projection  onto  a  two-dimensional  plane 
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along  a  predefined  direction. 

Another  application  where  intersections  occur  is  in  Constructive  Solid  Geometry 
(CSG)  where  an  object  is  constructed  from  sin:iple  building  blocks  by  means  of  inter- 
section and  union.  The  object  is  then  represented  by  the  tree  whose  leaves  are  the 
building  blocks  and  each  inner  node  stands  either  for  the  union  or  the  intersection  of 
the  objects  defmed  in  its  subtrees.  A  view  of  the  object  can  be  computed  by  post- 
order  traversal  of  the  defining  tree.  A  special  case  in  which  .\lgorithm  1  is  most  effec- 
tive is  when  the  object  is  simply  the  union  (or  intersection)  of  many  (polyhedral) 
building  blocks,  or  when  its  CSG  tree  has  at  most  two  levels.  For  an  arbitrarily 
defined  object,  however,  Algorithm  1  may  not  be  very  efficient. 


4.    Translating  a  Polyhedron  in  Three  Dimensions 

An  object,  B,  in  some  space  cannot  be  moved  to  any  arbitrary  position  if  there  are 
obstacles  present  which  it  has  to  avoid.  The  motion  planning  problem  for  B  is  to  cal- 
culate the  space  of  all  placements  of  B,  called  the  free  placements  of  B.  in  which  it 
does  not  collide  with  any  obstacle  (see  HSS^  for  a  recent  compendium  of  work  done 
on  motion  planning).  In  this  section  we  consider  special  cases  of  motion  planning  in 
which  B  is  allowed  to  translate  but  not  to  rotate.  The  problems  that  we  address 
make  sense  in  arbitrary  dimensions  but  for  simplicity  and  also  because  it  is  the  most 
important  setting,  we  discuss  only  the  three-dimensional  case.  The  much  simpler 
two-dimensional  case  has  been  studied  in  ;KLPS|,  [LS;,  ';PSSi,  JGSSJ.  If  the  object  as 
well  as  the  obstacles  are  polyhedra,  that  is,  their  boundaries  are  piecewise  linear,  then 
these  motion  planning  problems  lead  to  certain  envelope  questions  as  we  will  see 
below. 

Let  B  be  an  open  three-dimensional  polyhedron  bounded  by  k  facets  and  let 
Ai,A2,---,^m  b^  closed  convex  polyhedral  obstacles  bounded  by  a  total  number  of  n 
facets.  It  is  not  essential  that  B  is  open,  only  the  description  of  our  results  is  slightly 
easier  this  way  because,  otherwise,  we  have  to  allow  fi's  boundary  to  intersect  the 
boundary  of  an  obstacle  -  only  the  interiors  have  to  be  disjoint.  .\11  results  are  true 
for  B  closed  if  one  changes  the  phrasing  of  the  results  accordingly.  We  assume  that 
the  A,  are  convex;  so  any  non-convex  obstacle  is  split  into  convex  pieces  beforehand 
(see  C  ).  Our  goal  is  to  calculate  the  set  of  all  translates  of  B  that  avoid  the  obsta- 
cles. The  standard  approach  to  solving  this  problem,  initially  proposed  in  'LW',  uses 
Minkowski  differences  between  the  ,4,  and  B.  A  translate  B'  of  B  is  determined  by 
its  translation  vector  6,  that  is, 

B'  =  B+b  =  \x+b\x£B\. 

We  will  interchangeably  think  of  6  as  a  vector  and  as  a  point.  B'  intersects  obstacle 
A,  if  and  only  if  there  is  a  point  y£A,  and  a  point  x£B  such  that  x+b  =y,  which  is 
equivalent  to  b  =y—x.    Another  way  to  say  this  is  that  B'  and  A,  are  disjoint  as  long 
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as  b  is  not  in  the  set 


K,  =A,-B  =  \y-i\yeA,  and  xEB}. 


K,  is  known  as  the  Minkowski  difference  of  A,  and  B  and  is  sometimes  referred  to  as 
the  "expanded"  or  "grown"  obstacle.    It  is  clear  that  B'  lies  in  a  free  position  if  and 

n 

only  if  b  does  not  belong  to  /C  =  UA',.    We  can  thus  represent  the  set  of  free  posl- 

1  =  1 

tions  by  its  complement,  K' . 

To  get  a  handle  on  the  combinatorial  complexity  of  /C^  assume  that  obstacle  A^ 

m 

is  bounded  by  n,  facets;  thus,   S  n,  =n.    Except  for  degenerate  cases,  each  facet  of  /C, 

1  =  1 

is  the  Minkowski  difference  of  a  facet  of  A^  and  a  vertex  of  5,  of  an  edge  of  A,  and 
an  edge  of  fl,  or  of  a  vertex  of  A,  and  a  facet  of  B.  The  number  of  such  pairs  is 
Oik-n,),  in  contrast  to  the  planar  case  where  this  number  is  only  0(k  +  n,)  (see 
GRS  ).  This  implies  that  K,  is  bounded  by  at  most  0(A:n, )  faces.  As  a  matter  of 
fact,  the  number  of  facet-vertex  and  vertex-facet  pairs  is  0(A;  +  m, ),  only  the  number 
of  pairing  edges  can  be  quadratic.  Thus,  the  K,  altogether  can  be  modeled  by  O(k-n) 
triangles  which  implies  that  the  number  of  faces  needed  to  describe  K' ,  the  set  of  free 
placements  of  B,  is  0{k  n  ).  If  fl  is  non-convex  there  are  cases  where  the  boundary 
of  K'^  consists  of  Q{k  n  )  faces  (see  Figure  4.1:  the  ""triple  fork"  of  size  k  can  be 
locked  into  the  "three-sided  cage"  of  size  n  in  ri{k  n  )  different  ways)  which  shows 
that  the  trivial  bound  is  asymptotically  tight.  If  fl  is  convex,  then  no  such  example  is 
known  and  a  plausible  conjecture  is  that  K'  is  bounded  by  at  most  0{k^-n^-a(kn)) 
faces.  It  is  rather  easy  to  give  examples  where  the  complexity  in  question  is  Q{k^-n'). 
Our  goal  is  to  show  that  the  complexity  of  K'^  is  much  lower  than  proportional  to 
k  n  in  certain  important  cases,  or  failing  that,  to  show  that  the  complexity  of  a  sin- 
gle connected  component  of  K'  (which  is  often  all  we  need  to  consider)  is  small. 


Figure  4.1. 
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Consider  first  the  general  case.  If  we  are  interested  in  the  set  of  free  placements 
that  can  be  reached  by  B  from  its  initial  position  without  ever  interfering  with  obsta- 
cles, then  we  actually  ask  for  the  connected  component  of  K'^  that  contains  the  origin 
-  rather  than  for  the  entire  K'^ .  To  get  an  upper  bound  on  the  combinatorial  com- 
plexity of  this  connected  component  we  can  use  Theorem  4  in  [PS:  which  shows  that 
there  are  at  most  0{{k-n)^~^''^^)  faces  in  its  boundary.  Unfortunately,  the  proof  of 
this  result  is  non-constructive  and  does  not  lead  to  an  algorithm  that  constructs  the 
connected  component  in  o{k  -n  )  time.  Hence,  in  the  general  case,  even  though  this 
result  sheds  light  on  the  problem  structure,  no  satisfactory  solution  is  yet  available. 
Some  special  cases  of  this  general  problem  are  now  being  investigated  in  ,ASj. 

Things  are  much  improved,  however,  when  we  consider  the  special  case  in  which 
the  obstacles  A,  collectively  form  a  so-called  polyhedral  terrain.  This  is  a  piecewise 
linear  surface  that  intersects  every  vertical  line  in  exactly  one  point.  B  is  still 
assumed  to  be  an  arbitrary,  thus  not  necessarily  convex,  polyhedron  and  we  wish  to 
Hnd  all  free  placements  of  B  above  this  terrain,  £".  Again,  we  represent  such  a  place- 
ment B'  by  the  point  b  such  that  B'  =  B+b  and  use  the  preceding  analysis  to  obtain 
the  space  of  free  placements  A''^,  with  K  defined  as  above.  Of  course,  in  this  res- 
tricted case  the  resulting  set  /C^  is  connected:  B  can  be  translated  from  any  free 
placement  to  any  other  in  a  canonical  manner  by  first  moving  upwards  to  a  suffi- 
ciently high  position,  then  translating  horizontally  to  a  position  above  the  target  posi- 
tion and,  finally,  descending  to  the  desired  position.  Nevertheless,  the  calculation  of 
K'  is  significant  in  certain  applications.  Such  an  example  occurs  when  B  is  required 
to  maintain  a  fixed  maximal  vertical  distance  from  E,  for  example,  when  B  surveys  i7 
from  close  distance  as  it  flies  over  it.  Also,  by  preprocessing  the  boundary  of  K'  into 
a  data  structure  which  supports  fast  point  location  queries,  we  can  decide  in  loga- 
rithmic time  whether  or  not  a  given  placement  of  B  is  free  (see  [El,  chapter  IL  for  an 
optimal  data  structure  that  supports  point  location  queries).  This  method  can  also  be 
used  to  determine  the  point(s)  of  contact  of  B's  closure  with  S  as  it  is  lowered  until  it 
touches  E.  If  the  obstacles  together  form  a  polyhedral  domain,  then  the  boundary  of 
K'  is  the  upper  envelope  of  the  expanded  obstacles  K,  =A,—B.  Theorem  I'  in  [PSi 
and  the  algorithmic  results  in  Section  2  of  this  paper  now  imply  the  following  result. 

Theorem  4.1.  Let  5  be  a  polyhedron  bounded  by  k  facets  and  let  i7  be  a 
polyhedral  terrain  with  n  facets.  Then  the  number  of  faces  bounding  K' , 
the  set  of  free  placements  of  B,  is  0{k"-n''-a(k-n))  and  it  can  be  constructed 
in  the  same  amount  of  time. 

The  preceding  arguments  can  be  generalized  to  cases  where  we  allow  B  to 
translate  only  along  a  single  half-line  from  its  current  position.  Two-dimensional 
variants  of  this  problem  have  been  studied  extensively  by  Toussaint  and  others  (see 
[Tt|  for  a  survey).  By  applying  an  appropriate  projective  transformation  (similar  as 
in  Section  3)  we  can  assume  that  B  is  initially  at  infinity  and  is  allowed  to  descent 
along  vertical  lines  only.    For  an  arbitrary  collection  of  obstacles  A^  we  only  need  to 
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find  the  upper  envelope  of  the  expanded  obstacles  A",  =A,— fl  (again,  remennbering  to 
add  the  transformed  image  of  the  "'background"  plane  at  infmity).  For  each  vertical 
line,  this  envelope  gives  the  first  obstacle  hit  by  B  if  it  moves  along  this  line.  If  no 
obstacle  is  in  the  way  of  fl's  vertical  movement,  then  the  envelope  shows  the  former 
infinite  plane  as  the  obstacle  hit  first.  In  the  untransformed  space  this  corresponds  to 
B  moving  to  infinity  without  ever  hitting  an  obstacle.  Using  our  combinatorial  and 
computational  knowledge  about  envelopes  we  get  the  following  result. 

Theorem  4.2.  Let  fi  be  a  polyhedron  bounded  by  k  facets  and  let  Ai,A2,.-,Ar^ 
be  (possibly  intersecting)  convex  obstacles  bounded  by  a  total  of  n  facets. 
The  set  of  free  placements  of  B  that  are  reachable  by  translating  B  along  all 
possible  half-lines  is  bounded  by  0(a  Q(n))  faces  and  can  be  constructed  in 
0{n'a{n))  time. 

We  conclude  this  section  with  an  argument  that  supports  our  conjecture  that  the 
combinatorial  complexity  of  K'^  in  only  nearly  quadratic  in  kn  if  B  is  convex.  Here 
we  do  not  assume  that  the  obstacles  form  a  polyhedral  terrain.  We  show  that  this  is 
true  if  B  is  a  line  segment.  To  show  this  it  suffices  to  demonstrate  that  there  are 
only  O(n^)  free  placements  of  B  in  which  it  simultaneously  makes  contact  with  three 
obstacles.  These  triple  contacts  correspond  to  the  vertices  of  K' .  In  each  such  triple 
contact  one  of  the  contacts  must  be  at  a  point  of  B  different  from  its  two  endpoints. 
This  point  must  touch  an  edge  of  an  obstacle  provided  we  ignore  degenerate  cases. 
Fix  such  an  obstacle  edge  e  and  consider  the  plane  h  through  e  that  is  parallel  to  B. 
When  B  translates  within  plane  h,  maintaining  contact  with  e,  it  can  reach  at  most 
0(n)  placements  at  which  it  makes  two  more  contacts  with  the  obstacles  (see 
KLPS  ).    From  this  the  claim  follows  readily. 


5.    Stabbing  Line  Segments  and  Polytopes 

Finding  transversals  of  a  finite  set  of  objects  is  the  first  of  two  problems  discussed  in 
this  paper  that  relate  to  envelopes  by  means  of  a  geometric  transformation.  The 
second  such  problem  deals  with  certain  V'oronoi  diagrams  for  sets  of  point  clusters.  In 
Section  6  we  show  that  this  problem  is  in  fact  closely  related  to  the  stabbing  problem 
of  this  section. 

Let  S  be  a  finite  set  of  connected  objects  in  d  +  l  dimensions.  (We  use  d  +  l  to 
denote  the  dimensionality,  rather  than  d,  in  order  to  be  consistent  with  the  notation 
in  Section  2.)  A  hyperplane  is  a  transversal  of  S  if  it  intersects  each  object  in  5;  we 
also  say  that  it  stabs  S.  Since  a  hyperplane  intersects  a  connected  object  if  and  only 
if  it  intersects  its  convex  hull,  we  can  assume  without  loss  of  generality  that  all 
objects  in  S  are  convex.  We  consider  the  problem  of  finding  all  transversals  of  5,  or 
a  representation  of  this  set,  assuming  that  S  is  a  collection  of  convex  polytopes.    The 
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complexity  of  a  solution  will  be  measured  in  terms  of  n,  the  total  number  of  vertices 
of  the  polytopes.  In  three  dimensions,  Euler's  relation  implies  that  n  is  proportional 
to  the  number  of  edges  and  facets  bounding  the  polytopes.  This  is  no  longer  true  in 
four  or  higher  dimensions.  For  this  reason  we  restrict  most  of  our  discussion  to 
d+l  =3  dimensions  and  comment  on  the  difficulties  one  encounters  in  four  and  higher 
dimensions  at  the  end  of  this  section. 

Earlier  results  on  this  problem  can  be  found  in  E3!  which  gives  an  O(Mlogn) 
time  algorithm  for  5  a  set  of  n  line  segments  in  the  plane,  and  in  [AD  which  gives  an 
0(rt'')  algorithm  for  n  line  segments  in  <f+l=3  dimensions.  Indeed,  using  known 
algorithms  for  constructing  arrangements  of  hyperplanes  0(n  "^  )  time  algorithms  are 
possible  even  for  the  more  general  case  of  polytopes  in  ^  +  l>3  dimensions  (see  El, 
chapter  7  ).  We  show  in  this  section  that  this  straightforward  bound  can  be  improved 
to  0{n')  in  the  case  of  line  segments  and  to  0{n'a{n))  in  the  case  of  convex  polytopes 
in  three  dimensions.  These  results  are  optimal  in  the  worst  case  in  a  sense  that  will 
become  clear  later. 

Our  development  is  based  on  a  dual  transform,  P.  that  maps  a  point  to  a  plane 
and  vice  versa.  If  p  =(7''i,?i"2i^3)  is  a  point  in  three  dimensions,  then  we  define  the 
plane 

Notice  that  D(p)  is  non- vertical,  that  is,  intersects  the  i:3-axis  in  a  unique  point.  If  k 
is  a  non-vertical  plane  we  define  V{h)  =  p  such  that  h=D[p).  Thus.  D  is  involutary 
by  definition.  It  is  fairly  easy  to  show  that  D  preserves  incidence  relations  [pEh  if 
and  only  if  D{h)£D{p))  and  order  relations  (p  lies  vertically  above  h  if  and  only  if 
D{h)  lies  vertically  above  V(p)). 

We  next  extend  D  to  point  sets  and,  in  particular,  to  polytopes.  For  P  a  point  set 
in  three  dimensions  we  define 

D[P)  =  U  P(x), 

that  is,  D{P)  is  the  set  of  all  points  that  belong  to  at  least  one  plane  dual  to  a  point  of 
P.  We  call  D{P)  the  stabbing  region  of  P.  Since  D  preserves  incidences  we  have 
p^V(P)  if  and  only  if  plane  D{p)  inters£cts  P.  Figure  5.1  illust,rates  these  definitions. 
It  shows  P  as  a  convex  pentagon  in  two  dimensions  and  displays  the  stabbing  region 
of  P.  It  also  shows  a  line  intersecting  the  pentagon  and  it  dual  point  which,  of  course, 
belongs  to  D{P).  If  P  is  a  (connected)  polytope,  then  D{P)  is  the  set  of  all  points  that 
are  neither  vertically  below  all  planes  corresponding  to  vertices  of  P  nor  vertically 
above  all  such  planes.  This  is  because  for  such  a  point  x,  its  dual  plane,  P(x),  stabs  P 
and  thus  must  separate  at  least  one  pair  of  vertices  of  P.  For  P  a  line  segment  there 
are  only  two  vertices  and  thus  only  two  dual  planes.  Consequently,  D[P)  is  the  double 
wedge  of  points  that  lie  vertically  between  the  two  planes  corresponding  to  the  line 
segment's  two  endpoints. 
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Figure  5.1. 

By  definition,  a  transversal  is  a  plane  that  cuts  all  polytopes  in  5.  It  follows  that 
plane  h  is  a  transversal  of  5  if  and  only  if  its  dual  point,  D{h),  belongs  to 

s{s)  =  noiP). 

This  intersection  is  termed  the  stabbing  region  of  S.  It  is  a  representation  of  all 
transversals  of  S.  Notice  that  the  transformation,  as  currently  defined,  excludes  vert- 
ical planes  which  thus  have  to  be  treated  separately.  One  way  to  do  this  is  to  verti- 
cally project  the  polytopes  onto  the  plane  spanned  by  the  Xj-  and  Ji'^-xis  and  to 
solve  a  two-dimensional  stabbing  problem  for  this  set.  Every  transversal,  which  is 
now  a  line,  corresponds  to  a  vertical  transversal,  a  plane,  of  the  original  set,  S. 

For  a  given  polytope  P,  D{P)  is  the  set  of  points  below  or  on  Up  and  above  or  on 
{.?,  where  Up  (resp.  1?)  is  the  upper  (lower)  envelope  of  the  planes  dual  to  the  ver- 
tices of  P.  These  are  piecewise  linear  bivariate  functions.  Thus,  the  stabbing  region, 
${S),  is  the  set  of  points  between  or  on  the  lower  envelope  of  all  functions  Up  and  the 
upper  envelope  of  all  functions  Lp.  Using  results  from  Section  2  we  can  give  bounds 
on  the  combinatorial  complexity  of  S(S)  and  on  the  amount  of  time  needed  to  con- 
struct it. 

In  order  to  analyze  S(S)  we  model  each  function  Up  and  JLp  by  a  collection  of  tri- 
angles in  three  dimensions.  If  m  is  the  number  of  vertices  of  P,  then  Up  and  Lp  can 
be  decomposed  into  0{m)  triangles.  This  puts  us  into  the  situation  described  in 
Theorem  2.2  (i).  The  upper  bound  on  the  combinatorial  complexity  can  be  improved 
from  0(n^a(n))  to  O(a^)  if  5  is  a  set  of  n  line  segments.  This  is  because  each  Up  and 
Lp  is  composed  of  two  half-planes  that  meet  at  a  common  line.  The  improvement 
follows  now  from  Theorem  2.2  (ii). 

Theorem  5.1.    Let  5  be  a  set  of  convex  polytopes  in  three  dimensions  and  let  n 
be  the  total  number  of  vertices. 

(i)  The  number  of  faces  bounding  S{S)  is  0{n''a{n))  and  so  is  the 
amount  of  time  needed  to  construct  S{S).  Both  bounds  are  tight  in  the 
worst  case. 

(ii)    The  number  of  faces  bounding  S{S)  is  O(n^)  if  all  polytopes  in  5 
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are  line  segments.    In  this  case,  0(n  )  time  suffices  to  construct  S{S).    Both 
bounds  are  tight  in  the  worst  case. 

Using  the  lower  bounci  examples  indicated  in  'PS]  it  is  not  difficult  to  prove  that  all 
bounds  are  asymptotically  tight  in  the  worst  case.  In  this  context  it  is  interesting  to 
note  that  examples  with  Q{n  ci{n))  faces  can  be  modeled  even  with  the  restriction 
that  all  polytopes  in  5  are  triangles  in  three  dimensions. 

If  we  specialize  the  computational  results  of  Theorem  5.1  to  d  +  l  =2  dimensions, 
we  get  an  algorithm  that  constructs  the  stabbing  region  of  polygons  with  a  total  of  n 
vertices  in  0{na{n)\ogn)  time,  and  in  0{nlogn)  time  if  all  polygons  are  line  segments. 
The  former  result  is  new  although  it  follows  easily  from  the  combinatorial  analysis  of 
two-dimensional  envelopes  in  [HS]  and  the  algorithmic  techniques  in  iE3;.  The  latter 
result  dates  back  to  E3i.  Note  the  logn  term  in  the  time-complexity  that  comes  up 
in  two  dimensions.  The  reason  for  this  extra  term  is  that  the  homogeneous  solution 
of  the  recurrence  relation  that  describes  the  time-complexity  is  essentially  the  same  as 
the  additive  term  (see  Section  2).  In  three  dimensions  the  additive  term  is  signifi- 
cantly larger  than  the  homogeneous  solution  which  explains  why  the  log  a  term  disap- 
pears. 

We  conclude  this  section  with  a  few  remarks  about  the  generalization  of  our 
methods  to  d  +  l>A  dimensions.  The  first  difficulty  that  arises  is  combinatorial  and 
concerns  the  decomposition  of  the  Up  and  Lt>  into  <f-simplices.  The  total  number  of 
vertices  of  the  input  polytopes,  n,  is  proportional  to  the  number  of  facets  of  the  Up 
and  {.  ?,  but  it  might  very  well  be  that  the  number  of  lower-dimensional  faces  of  the 
Up  and  Lp  by  far  exceeds  0(n).  Indeed,  in  <f+l=4  dimensions  their  number  is  n(a') 
if  S  consists  of  a  constant  number  of  cyclic  polytopes.  Such  constellations  need  more 
than  0(n)  simplices  to  model  the  boundary  of  the  stabbing  region  by  two  envelopes  of 
simplices.  This  weakens  our  bounds  on  the  combinatorial  complexity  of  stabbing 
regions  which  use  envelope  bounds.  Even  if  we  had  a  method  that  circumvents  the 
sketched  difficulty,  there  is  no  algorithm  known  that  constructs  the  stabbing  region  in 
time  o(rt  *  )  because  of  reasons  explained  in  Section  2.  But  0{n  )  is  straightfor- 
ward if  one  uses  arrangement  algorithms  as  mentioned  above. 


8.    Voronoi  Diagrams  of  Point  Clusters 

In  this  section  we  consider  applications  of  envelopes  to  a  certain  generalization  of 
Voronoi  diagrams.  This  generalization  can  be  defined  in  any  number  of  dimensions, 
and  we  will  do  so,  but  our  discussion  of  the  combinatorial  and  computational  com- 
plexity is  mostly  confined  to  the  plane.  The  diagram  that  we  have  in  mind  bears 
close  relationship  to  the  notion  of  complete  linkage  clustering  (see  e.g.  Hi).  For  this 
clustering  method,  the  distance  between  any  two  clusters  is  defined  as  the  maximum 
distance  between  any   two  points,  one  of  each  cluster.     We  will  come  back  to  this 
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clustering  method  at  the  end  of  this  section. 

Let  6  denote  the  Euclidean  distance  function.  For  a  set  of  points,  C,  called  a 
cluster,  and  for  a  point  p,  define 

8(p,C)  =  mz.x\8{p,x)\zeC] 

as  the  distance  between  p  and  C.  In  most  cases  we  let  C  be  finite  and,  whenever  it  is 
convenient  to  have  C  infinite,  it  will  be  the  convex  hull  of  a  finite  number  of  points 
in  which  case  the  maximum  of  the  distances  between  p  and  points  of  C  is  well 
defined.  The  definition  of  6{p,C)  implies  that  the  closed  ball  with  center  p  and 
radius  S{p,C)  contains  cluster  C.  In  fact,  it  is  the  smallest  ball  centered  at  p  for 
which  this  is  true.  Let  now  S  be  a  finite  set  of  clusters.  The  Voronoi  diagram  of  5, 
V{S),  is  a  partition  of  space  into  maximal  (but  not  necessarily  connected)  domains, 
one  for  each  cluster,  such  that  a  point  p  belongs  to  the  domain  associated  with  cluster 
C  if  6{p,C)<6{p,D)  for  all  clusters  D  ^C  in  S.  See  Figure  6.1  for  an  example.  It  is 
convenient  to  replace  a  cluster  C  by  its  convex  hull  which  is  all  right  since  the  dis- 
tance of  a  point  p  from  C  is  the  same  as  its  distance  from  the  convex  hull  of  C .  The 
clusters  in  Figure  6.1  are  A  (two  points),  B  (two  points),  and  C  (three  points).  The 
domain  of  B  consists  of  two  connected  components  which  we  call  regions.  The  solid 
lines  show  the  Voronoi  diagram  and  the  dashed  lines  decompose  each  region  using  the 
furthest  point  Voronoi  diagram  of  the  cluster.  This  is  the  diagram  that  associates 
with  each  point  the  part  of  the  domain  for  which  the  point  is  the  furthest  point  of  the 
cluster.  The  significance  of  this  decomposition  is  that  it  shows  which  point  of  the 
cluster  attains  the  distance  to  the  cluster  and  where  it  does  so. 

Below,  we  discuss  some  properties  of  this  kind  of  Voronoi  diagram.  First,  we 
demonstrate  that  they  are  closely  related  to  envelopes  and  thus  derive  general  upper 
bounds  on  their  complexity.  Second,  we  study  the  special  case  where  the  convex  hulls 
of  the  clusters  are  disjoint.  It  turns  out  that  this  condition  reduces  the  combinatorial 
complexity  of  the  diagrams  dramatically.  For  simplicity,  we  restrict  our  attention  to 
the  two-dimensional  case.  Three-  and  higher-dimensional  cases  can  be  treated  in  the 
same  way  as  the  stabbing  problem  in  four  and  higher  dimensions  (see  Section  .5);  we 
thus  omit  all  details  pertaining  to  these  extensions.  Specifically,  we  prove  the  follow- 
ing theorem. 

Theorem  8.1.  Let  S  be  a  set  of  clusters  in  the  plane  and  let  n  be  the  sum  of 
the  cardinalities  of  the  clusters. 

(i)   The  number  of  faces  of  'V(S)  is  0(n^a(n)). 

(ii)  If  each  cluster  consists  of  one  or  two  points,  then  the  number  of 
faces  of  'V(S)  is  0{n^),  and  this  is  tight  in  the  worst  case. 

(iii)  If  the  convex  hulls  of  any  two  clusters  are  disjoint,  then  V(5)  con- 
tains at  most  \S\  regions.  The  number  of  edges  and  vertices  in  this  case  is 
0(n). 


20 


B's  domain 


.■^     A'8    domai 


C's  domain 


Figure  8.1. 

Proof.  The  proof  consists  of  three  fairly  independent  steps.  First,  we  demonstrate 
the  upper  bounds  in  (i)  and  (ii)  by  means  of  envelopes  in  three  dimensions.  Second, 
we  construct  an  example  that  proves  the  lower  bound  on  the  maximum  complexity 
stated  in  (ii).  Finally,  we  show  that  the  domain  of  a  cluster  is  connected  if  the  convex 
hulls  of  any  two  clusters  are  disjoint.    This  will  lead  to  a  proof  of  (iii). 

In  order  to  relate  V(S)  to  an  envelope  of  triangles  in  three  dimensions  we  use  two 
geometric  transforms.  Let  p  =(7r,,7r2)  be  a  point  in  the  two-dimensional  plane  which 
we  associate  with  the  plane  13=0  in  three  dimensions.  The  Hrst  transform,  U ,  pro- 
jects  p  vertically  onto  the  paraboloid  of  revolution  given  by  Xj^xf+Xj  ,  that  is. 

The  second  transform,  £,  maps  p  to  the  unique  plane  that  touches  the  paraboloid  in 
point  ^[j>\  that  is, 

£(p):  X3=27riXi+27r.,J2-('^f+^2)- 
These  transforms  can  be  used  to  express  distance  information  in  two  dimensions  as 
combinatorial  information  in  three  dimensions.  The  crucial  property  here  is  that 
(5^(p,x),  the  square  of  the  distance  between  points  p  and  x  in  the  plane,  is  equal  to  the 
vertical  distance  between  point  0'(x)  and  the  vertical  projection  of  x  onto  plane  f(p) 
(see  ;E1,  chapter  l]).    The  distance  from  x  to  a  point  p  is  thus  the  square  root  of  the 
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vertical  distance  from  point  C'(x)  down  to  plane  S{p).  Consider  now  a  cluster  of 
points,  C,  instead  of  a  single  point.  The  distance  from  z  to  C,  6(x,C),  is  equal  to  the 
square  root  of  the  largest  vertical  distance  from  U{x)  to  any  of  the  planes  £{p),  p£C, 
that  is,  to  the  lower  envelope,  Lc,  of  all  these  planes.  If  C  contains  m  points,  then 
this  envelope  is  the  boundary  of  a  convex  polyhedron  consisting  of  at  most  m  facets 
which  can  be  decomposed  into  0(m)  triangles.  Let  .V/  be  the  upper  envelope  of  all 
the  surfaces  Lq^  C£S.  Then  a  point  i  lies  in  the  domain  of  cluster  C  exactly  when 
Lq  is  vertically  nearest  to  U{x)  among  all  Lq,  DG5.  Since  U{x)  lies  above  all  sur- 
faces Ld,  this  is  equivalent  to  Lq  attaining  the  upper  envelope  M  at  r.  Thus,  we  get 
^(S)  by  constructing  the  upper  envelope  of  the  surfaces  Lq,  or  more  specifically  of 
the  triangles  composing  these  surfaces,  and  then  projecting  the  faces  of  the  envelope 
vertically  onto  the  plane  X3=0.  The  upper  bound  of  (i)  follows  now  immediately 
from  Theorem  l'  in  ;PS'.  If  each  cluster  consists  of  at  most  two  points  each,  the  Lp 
are  either  single  planes  or  two  half-planes  glued  together  along  a  common  line.  For 
these  functions  we  have  an  upper  bound  of  O(n')  for  the  combinatorial  complexity  as 
argued  in  Section  2.    This  proves  the  upper  bound  in  (ii). 

We  next  show  the  lower  bound  in  (ii).  Assume  without  loss  of  generality  that 
n  =4k  for  some  integer  k.  We  describe  a  set  S  of  2k  clusters,  each  containing  two 
points,  such  that  ^S)  has  Q{n  )  regions.  This  example  is  illustrated  in  Figure  6.2.  It 
consists  of  k  point  pairs  on  the  j:[-axis  such  that  the  i  pair  can  be  obtained  by  mov- 
ing the  (i— 1)^'  pair  a  distance  t  to  the  right,  for  2<i  <A:  and  e>0  sufficiently  small. 
More  specifically,  we  choose  the  first  pair  of  points  at  locations  (—1,0)  and 
(1— (A:  — l)-e,0),  with  ke<l/2.  Note  that  this  implies  that  the  X2-axis  is  a  symmetry 
axis  of  the  2k  points.  The  Voronoi  diagram  of  these  clusters  only  consists  of  A:— 2 
vertical  slabs  of  width  i,  one  for  each  pair  except  for  the  first  and  the  last,  and  two 
half-planes,  one  for  the  first  point  pair  and  one  for  the  last.    .\dd  how  another  k 


Figure  6.2. 
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point  pairs  on  the  2;.,-axis.  To  describe  how  these  points  are  to  be  chosen,  we  take  a 
point  p  and  move  it  upwards  inside  one  of  the  slabs.  For  each  location  of  p  we  con- 
sider the  smallest  disk  with  center  p  that  contains  the  horizontal  point  pair 
corresponding  to  the  slab.  By  construction,  the  disk  contains  no  other  horizontal 
cluster.  Since  the  horizontal  clusters  are  almost  identical,  the  wandering  disk  is 
almost  the  same  for  point  p  moving  in  any  other  slab.  Thus,  we  can  choose  point 
pairs  on  the  j2-axis  such  that  the  disk  alternates  between  containing  one  vertical 
cluster  and  containing  no  vertical  cluster.  Each  slab  is  then  decomposed  into  2k +1 
pieces  k  +  l  of  which  defme  the  domain  of  the  corresponding  horizontal  cluster.  The 
lower  bound  in  (ii)  follows. 

Finally,  we  prove  that  each  domain  is  connected  if  the  convex  hulls  of  the  clus- 
ters are  pairwise  disjoint.  The  proof  takes  two  steps.  First  it  verifies  that  the  skele- 
ton of  the  furthest  point  Voronoi  diagram  of  a  cluster,  for  short  the  skeleton  of  the 
cluster,  intersects  the  domain  of  the  cluster  in  a  connected  tree.  Second  it  shows  that 
an  arbitrary  point  of  the  cluster  can  be  connected  to  this  tree  by  a  straight  line  seg- 
ment that  lies  entirely  within  the  domain  of  the  cluster.  We  do  the  second  step  first 
because  it  is  simpler  than  the  first  step.  Let  p  be  a  point  that  belongs  to  the  domain 
of  a  cluster  C  (see  Figure  6.1).  By  definition,  C  is  the  only  cluster  that  is  fully  con- 
tained in  the  closed  disk  with  center  p  and  radius  6{p,C).  Let  c  be  the  point  in  C 
such  that  6{p,c)  =  6{p,C).  If  c  is  not  unique,  then  p  already  belongs  to  the  skeleton 
of  C.  Otherwise,  move  p  straight  towards  c  until  it  runs  into  an  edge  of  the  skeleton. 
The  disk  at  every  intermediate  location  of  p  lies  strictly  inside  p's  original  disk  which 
implies  that  it  contains  no  other  cluster  besides  C.  It  is  clear  that  p  must  run  into  an 
edge  of  the  skeleton  for,  otherwise,  the  disk  of  p  would  eventually  vanish,  which  can 
only  mean  that  c  is  the  only  point  in  C  and  thus  is  equal  to  the  skeleton  of  C  by 
definition. 

We  now  prove  that  if  two  points,  x  and  y,  on  the  skeleton  of  C  belong  to  the 
domain  of  C,  then  all  points  between  x  and  y  also  lie  in  this  domain.  Note  that  the 
set  of  points  between  j  and  y  is  well  defined  since  the  skeleton  is  a  tree.  We  show 
this  in  two  steps.  Assume  first  that  x  and  y  lie  on  the  same  edge  of  the  skeleton  (see 
Figure  6.3(a)).  Let  z  be  a  point  between  x  and  y  and  let  d^,  d  ,  and  d,  be  the  smal- 
lest disks  around  i,  y  and  z  that  contain  C .  By  construction,  the  bounding  circles  of 
the  three  disks  have  two  common  points  which  are  the  two  points  in  C  that  maximize 
the  distance  to  z,  y,  and  z.  Thus,  the  line  segment  connecting  these  two  points,  call 
it  5,  is  a  common  chord  of  all  three  circles  and  also  belongs  to  the  convex  hull  of  C . 
Consequently,  d^  is  contained  in  the  union  of  d^  and  d  and,  by  assumption,  neither 
d^   nor  d    contain  a  cluster  besides   C.    Thus,  the  only  possibility  for  d,  to  contain 


'The  skeleton  of  a  cluster  C  is  the  set  of  points  i  such  that  there  are  at  least  two  points  c^C  that 
maximize  6{i,C).  It  is  a  straight  line  tree  with  at  most  m—2  vertices  and  2m-3  edges,  if  m  Is  the 
number  of  points  in  C.    For  convenience  we  define  the  skeleton  to  be  the  point  c  itself  if  C=\c\. 


23 


another  cluster  D^C  is  that  at  least  one  point  of  D  belongs  to  d^—d^C.d^  and  at 
least  one  other  point  belongs  to  d^-d^Qd^.  However,  any  two  such  points  are 
separated  from  each  other  by  5  (see  Figure  6.3(a))  contradicting  our  assumption  that 
the  convex  hulls  of  C  and  D  are  disjoint. 

In  the  second  step  we  show  that  if  x  and  y  lie  on  adjacent  edges  of  the  skeleton, 
then  the  common  vertex  of  these  two  edges  also  belongs  to  the  domain  of  C.  Call 
this  vertex  z  (see  Figure  6.3(b)).  The  proof  is  exactly  the  same  as  before.  We  still 
have  d^Qd^Udy  (as  will  be  shown  below)  and  the  regions  d^—d^  and  d^—dy  are  again 
separated  inside  d^  by  a  chord  (actually  two  chords)  connecting  two  points  of  C .  The 
line  segment  connecting  the  two  points  of  C  furthest  from  x  or  the  one  connecting 
the  two  points  furthest  from  y  are  such  chords.  Thus,  if  d^  contains  a  cluster  D,  then 
the  convex  hull  of  D  must  intersect  the  convex  hull  of  C . 

We  still  have  to  give  a  proof  that  d^  is  indeed  contained  in  the  union  of  d^  and 
dy.  Notice  that  the  bounding  circles  of  the  three  disks  all  go  through  a  common  point 
of  C  which  we  call  p.  (This  is  the  point  of  C  furthest  from  all  points  of  the  two  line 
segments  iz  and  yz.)  Consider  the  two  half-planes  that  contain  p  and  are  bounded 
by  the  lines  supporting  iz  and  yz.  By  construction,  neither  xz  nor  yz  lies  in  the 
boundary  of  the  intersection  of  those  two  half-planes.  Let  now  z'  be  the  point  that  is 
collinear  with  x  and  y  as  well  as  with  p  and  z  (see  Figure  6.3(b)).  Note  that  z'  lies 
between  i  and  y;  otherwise,  either  xz  or  yz  would  be  in  the  boundary  of  the  intersec- 
tion of  the  two  half-planes  defined  above.  By  the  same  reason,  z  and  p  lie  on  the 
same  side  of  the  line  through  x  and  y.  The  disk,  d^-,  bounded  by  the  circle  centered 
at  z'  that  goes  through  p  is  contained  in  d^ljdy.  Now,  d^C.d^<  unless  z  and  p  lie  on 
different  sides  of  z'.  But  this  is  impossible  because  z  lies  on  the  same  side  of  the  line 
through  X  and  y  as  point  p. 


(a) 


(b) 


Figure  6.3. 
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Thus  we  have  shown  that  the  domain  of  each  cluster  in  S  is  connected  (it  can  be 
empty,  though)  if  their  convex  hulls  are  pairwise  disjoint.  The  first  part  of  claim  (iii) 
follows.  This  does  not  imply  that  the  number  of  edges  and  vertices  of  l'(S)  is  also 
0(|S|).  By  Euler's  relation  for  planar  graphs  this  is,  however,  true  for  vertices  of 
V(5)  that  are  incident  to  at  least  three  edges.  .\  vertex  of  degree  two  must  also  lie  on 
an  edge  of  the  skeleton  of  the  cluster  in  one  of  the  two  adjacent  domains.  Since  each 
edge  of  the  skeleton  meets  the  boundary  of  its  cluster's  domain  in  at  most  two  points 
(this  follows  from  the  above  argument)  we  conclude  a  linear  upper  bound  on  the 
number  of  degree  two  vertices  and  thus  of  edges  of  V(S).  D 

We  remark  that  V(5)  can  be  further  refined  by  decomposing  each  region  of  ^(5) 
into  subregions  by  the  corresponding  skeleton.  This  yields  a  refined  convex  subdivi- 
sion of  the  plane  so  that  a  point  x  belongs  to  a  subregion  corresponding  to  point  c  of 
cluster  C  if  and  only  if  x  is  nearest  to  C  and  its  distance  to  C  is  attained  by  c.  In 
the  special  case  where  the  convex  hulls  of  the  clusters  are  disjoint,  the  domain  of  clus- 
ter C  is  thus  decomposed  into  at  most  \C\  subregions.  Euler's  relation  implies  that 
the  numbers  of  edges  and  vertices  of  the  refined  subdivision  are  0(n),  n  the  total 
number  of  points  in  all  clusters. 

To  construct  V(S)  for  an  arbitrary  collection  5  of  clusters  we  can  use  the  three- 
dimensional  envelope  algorithm  of  Section  2  which  implies  the  following  result. 

Theorem  8.2.  Let  5  be  a  set  of  clusters  where  n  is  the  sum  of  the  cardinalities 
of  the  clusters. 

(i)    V(5)  can  be  constructed  in  0{n  ci{n))  time. 

(ii)  If  each  cluster  is  of  size  one  or  two,  then  0(n")  time  suffices  and 
this  is  optimal  in  the  worst  case. 

There  are  two  major  open  problems  concerning  cluster  Voronoi  diagrams  that 
remain.  Is  0{n'a{n))  for  the  combinatorial  complexity  of  V(S)  tight?  A  better  upper 
bound  (maybe  0{n  ))  would  also  improve  the  time  bound  in  Theorem  6. 2(1).  Second, 
can  l'(S)  be  constructed  in  less  than  quadratic  time  (maybe  O(nlogfi))  if  the  convex 
hulls  of  the  clusters  are  pairwise  disjoint?  An  affirmative  answer  to  the  second  ques- 
tion could  also  be  relevant  to  complete  linkage  clustering  of  n  points  in  the  plane.  In 
this  method,  the  points  are  considered  to  be  individual  clusters  initially,  and  at  each 
stage  the  two  nearest  clusters  are  merged  until  all  points  belong  to  the  same  cluster. 
.\s  mentioned  at  the  beginning  of  this  section,  the  distance  between  two  clusters  is 
defined  as  the  maximum  distance  between  any  two  points,  one  from  each  cluster. 
The  most  efficient  algorithm  known  for  this  problem  takes  0{n^)  time  and  0(n) 
storage  for  the  entire  sequence  of  merges  (see  iDf!).  At  each  stage,  the  set  of  points  is 
partitioned  into  a  collection  of  clusters.  The  Voronoi  diagram  of  these  clusters  is 
potentially  useful  since  the  two  nearest  clusters  also  have  a  common  edge  in  the 
diagram.  The  problem  is  now  to  maintain  the  diagram  through  a  sequence  of  n—l 
cluster  merges.    In  the  case  of  single  linkage  clustering,  where  the  distance  between 
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two  clusters  is  the  minimum  distance  between  any  two  points  one  from  each  cluster,  a 
similar  approach  yields  an  O(nlogM)  time  algorithm  (see  El  ).  The  latter  clustering 
method  is  intimately  related  to  the  notion  of  the  minimum  spanning  tree  of  the 
points. 
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